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Abstract: Reentry vehicles are required to meet stringent boundary conditions and path 
constraints, while minimizing certain cost functional. Parameter optimal control is easier and 
faster to solve than indirect methods for complex problems, such as the problem of re-entry 
trajectory optimization subject to multiple path constraints. In this paper, a direct transcription 
method is applied for finding optimal trajectory of 2 problems. Firstly, Newtonian method is 
used for finding aerodynamic coefficients of blunt nose conical re-entry capsule in hypersonic 
regime. Using these coefficients, trajectory satisfying path constraints corresponding to heat 
influx, G-load and control angles is generated while minimizing the vehicle manoeuvres. In 
the second problem Apollo capsule’s multiphase optimization problem is solved based on the 
aerodynamic coefficients obtained from wind tunnel test data. 
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1 . INTRODUCTION 

Indian Space Research Organization (ISRO) has success¬ 
fully carried out Space Reentry Experiment (SRE-1), 
where a blunt nosed conical capsule was subjected to reen¬ 
try and was recovered. SRE-1 orbited around the Earth 
in polar sun-synchronous orbit. It was de-boosted and 
recovered by landing in Bay of Bengal. Though this vehicle 
did not have trajectory control capability, future mission 
might well have it. This paper is oriented to contribute to 
such future projects. 

Reentry initial conditions can deviate from the nominal, 
which need to be compensated for reaching the final 
destination accurately. In such conditions, trajectory needs 
to be designed for meeting the terminal conditions and 
safety requirements, while minimizing fuel consumption 
due to vehicle manoeuvres. This problem has been solved 
in this work for a body of shape similar to SRE-1. 

This paper addresses 2 problems of reentry trajectory op¬ 
timization using General Purpose Optimization Software 
(GPOPS-II). The first one is single phase optimization of 
SRE capsule’s trajectory subject to path constraints. The 
second problem deals with a multiphase problem applied 
to reentry of Apollo capsule. 

2 . LITERATURE SURVEY 

The theory of direct transcription method has been ex¬ 
plained in many texts and the method is widely used re¬ 
cently for through numerous commercial and open source 
software. The similar problem as addressed in this paper of 
reusable launch vehicle has been solved in using Sparse 
Optimal Gontrol Software (SOGS) and is also solved using 


GPOPS in This paper presents a different version of 
the reentry problem and also shows a solution method of 
a more practical multiphase problem in which the various 
regimes are captured through the aerodynamic coefficient 
variations during vehicle journey. 

3 . RE-ENTRY DYNAMIGS 


This section deals with the governing equations of transla¬ 
tional motion of re-entering body, and the path constraints 
for trajectory design 


3.1 Equations of motion 


Eollowing are the translational equations of motion for 
reentry body in geographic frame for round rotating earth. 

0 Longitude 

(j) Latitude 

r Radius from earth centre 

V Velocity 

X Heading angle with respect to local north 

7 Elight path angle (above local horizontal) 

uje Earth’s angular speed 


If summation of all the forces are represented in wind 
frame as [Fi, F2, Fs], the 6 state equation of vehicle mass 
centre are: 


^ = ('T’/r) cos7sinx(cos(5) ^ 

(j) = {v/r) cos 7 cos X 
r = sin 7 

pp 

V = L0pcosS(— cos 7 cos X cos (5 + sinxcosJ) + — 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 
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X = 2co’e(— tan 7 COS x cos S + sin^) + ... 

• c- • c- 

-Sin X cos 0 sin 0 + ... 

V cos 7 

- COS 7 sin X tan o + -- 

r Mv cos 7 


7 = 2co’e sin X cos ^ H-^ cos^fsinxcosxsin^ + ... 

_L ^ I ^3 

COS 7 COS 0) H— cos 7 + —— 

' ^ r ' Mv 

3.2 Path constraints 


The re-entry vehicle over its journey through the atmo¬ 
sphere faces 2 major loads: Thermal heat flux and G- 
forces. 


Heat Load: The convective heat transfer to the body 
forms the major source of heat influx. A simplified rela¬ 
tionship for the stagnation point convective heat rate is 
obtained as : 

/ \ 0.5 

q = 1.45 X 10“^ ^ wattlcm^ (7) 

The limits for this loading are 60 Wjew? for a manned 
vehicle and 500 Wjew? for an unmanned vehicle. 



Fig. 1. Cone segmentation of axisymmetric body 
^.1 Aerodynamic coefficients for blunt cone 

If Newtonian method is applied to a sharp nose cone, 
the aerodynamic coefficients can be found out analytically. 
Following are their expressions in terms of semicone angle 
(Oc) and angle of attache (a)- 
Axial force coefficient: 

Ca = 2sm^0c + sin^(a(l — 3sin^6>c) (10) 


G forces: In absence of thrust, the only contributors to 
the body forces is the aerodynamics. The value of upper 
limits of these forces are 4g for manned vehicle and 20g 
for unmanned vehicle. 


< Frr 


Normal force coefficient: 

Cn = cos^ Oc sin^ 2(a 

Pitching moment coefficient: 

Cm = —2/3 (Cat/ tan 6>c cos^ 6>c ) 


( 11 ) 

( 12 ) 


Also, a reentry capsule owing to heat flux limitations, 
needs to orient, so as to have its blunt nose windward. This 
requirement imposes limitation on angle of attack. These 
limitations are imposed as the constraints when optimal 
trajectory is designed. 

4. HYPERSONIC AERODYNAMICS 

The hypersonic regime forms the longest range traversal 
phase during re-entry. Hence, optimization of trajectory 
during this phase shall be of high significance. In this 
section the method of finding aerodynamic coefficients for 
a capsule with blunt nose cone have been presented . 

According to Newtonian theory, a stream of particles move 
in a straight parallel path toward the surface. There is no 
random motion, usually associated with a microscopic gas 
particles. A fluid particle impinging on the surface of a 
body would lose all the momentum normal to the surface, 
whereas the momentum tangential to the surface would 
be retained. The force acting on the body is then found 
out by resultant rate of change of momentum transfer. 
Moreover, the leeward surfaces would not be acted upon 
by any forces, as they do not interact with flow. 

This theory matches closely with the flow in hypersonic 
regime, due to shock phenomenon. Eollowing are the 
pressure coefficients’ relation for flat surface: 

Cp = 2cos2 6 », C'/ = 0 (9) 

where, 0 = angle between local normal and wind velocity 
vector 


This method can be applied to any axisymmetric body, 
which can be considered to be made up of large number of 
conical segments, as shown in figure I. The coefficients for 
the entire body are obtained as the combination of those 
for individual cone segments as follows: 

N 

Ca = Cai X r^/rjf + ^[<^^*(4-4-1)/^^] (13) 

i=2 

N 

Cn = Cjvi X rl/rh + YiF^Mi - Pi)/^ n] (14) 

i=2 

N 

Cm = Cmi X rl/r% + -Pi)/P •••( 15 ) 

- CN,Ai - Piilsi - Xi)/r%] 

This method is applied to the reentry capsule of shape 
similar to SRE-I, as shown in figure 2. These coefficients 
are in body frame, which are transformed to aerodynamic 
frame to obtain lift, drag and pitching moment coefficients. 
The values are computed at various angles of attack from 
-30 to 30 degree. The optimal control method requires an 
analytical expression of these coefficient against aerody¬ 
namic angles. Hence, polynomial curve fitting was done 
using Matlab® to find following coefficients as a function 
of angle of attack. 

Cd = -0.9596 + 1.294 (16) 

Cl = 0.4451 - 0.6874 a (17) 

Cm = 0.4873 q;3 + 0.7377 _ 0.77 a-0.3861 (18) 
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aloha 

Fig. 3. Cl, Cd curve fits 

The variation in Cl, Cd is shown in figure 3. It is also 
observed that due to offset eg, the vehicle trims at angle 
of attack {atrim) = —21.6 degree. 

5. TRAJECTORY OPTIMIZATION 

In this section, first a general continuous time formulation 
of optimal control problem is presented, which is followed 
by the direct transcription method to convert the problem 
into a Non-Linear Programming problem (NLP), using 
Legendre-Gauss-Radau (LGR) pseudospectral method. 

5.1 Optimal control problem formulation 

For, state variable x{t) Initial time ri eR, Final time 
Tat 6 R, to determine control input vector u{r) e to 
minimize the objective 

J = <I>(x(ro),to,x(TAr),t/) + ... 

1 

j g[xiT),u{T),to,tf]dt 

-1 

Subject to: 


• Plant dynamics 

^ u{t), t; to, tf) e i?” (20) 

• Boundary conditions 

^{x{Ti),to,x{TN),tf) =0 e (21) 

• Path constraints 

C(x(t), i4(r), r; to, t/) < 0 e (22) 

5.2 Direct transcription formulation 

Solution of the above problem by analytical or indirect 
methods is infeasible for complex problems such as the 
one of trajectory optimization. Hence, the problem is tran¬ 
scribed to NLP, solving which is easier. The transcription 
is done at LGR points. Following are the transcribed 
variables of the original problem. 

• States: 

n 

(23) 

i=l 

where, Li{r) is the Lagrange polynomial. 

• State derivatives 


n 


2 — X 

(24) 

• Boundary conditions 


Y7v,TAr) = 0 

(25) 

• Path constraints 


C{Xk,Uk,rk;to,tf) < 0 

(26) 

• Dynamic constraints 


f2DkiXi-^-^^f{Xk,Uk,Tk;to,tf) 

2=1 

Y7 

o 

II 




where, D^i = Li{ri) is {N—1) xN Radau psuedospec- 
tral differentiation matrix. 

• Gost functional 

J = $(Xi, to, YAr,t/) + 

tf—to fV TT -I- -I- \ 

—-— 2^Wkg[Xk,Uk,rk;to,tf) 

k=l 

Equations 23 to 27 alongwith cost functional 28 form an 
NLP, whose solution is the solution of original continuous 
time optimal control problem. 

6. PROBLEM 1 DEFINITION AND RESULTS 

The SRE-1 capsule after de-boosting operation enters 
atmosphere at 100km altitude at 8 km/sec. The end of 
hypersonic phase occurs at about 36 Km, and is considered 
to take place over Bay of Bengal at 16 degree latitude, 85 
degree longitude. Considering initial flight path angle of 
-1 degree, the remaining nominal entry point states are 
obtained as follows by knowledge of terminal states and 
original orbit. 

Nominal initial conditions: 

Latitude = —18.85 deg. Longitude = 91.8 deg 


70 



























2014ACODS 

March 13-15, 2014. Kanpur, India 



Fig. 4. Altitude vs time 

Altitude = 100 km^ Speed = 8 km/s 

Heading = —11.19 deg, Flight path angle = —Ideg 

But it is possible that these nominal initial conditions are 
not achieved due to various reasons, and still the capsule 
is required to precisely reach the desired destination. The 
optimal control problem, thus is to find the trajectory 
such that with feasible deviations in initial conditions, the 
terminal point is reached with minimum fuel consumption. 
The deviations are considered in the longitude (1 degree), 
flight path angle (5%), and heading angle (5%). Since 
we are only considering translational dynamics, attitude 
dynamics and reaction control system is not modelled, and 
hence fuel consumption is not directly available. But, a 
measure of this is the deviation of angle of attack from 
trim angle of attack, which is considered as the optimality 
criterion. Hence, optimal control problem definition is to 
find angle of attack {a) and bank angle (/3) for state 
dynamics 1 to 6: 

Vehicle Properties'^]: 

Mass = 550 kg^ Area = 0.7854 sq.m 

Initial conditions: 

Latitude = —18.85 deg. Longitude = 92.8 deg 

Altitude = 100 km^ Speed = 8 km/s 

Heading = —11.75 deg, Flightpathangle = —1.05 deg 

Constraints: 

Angle of attaek = —30 deg to 30 deg 
Bank angle = —90 deg to 90 deg 
Max.heat rate = 60 WjevnS., 

Max. G load = Ag 

Terminal condition: 

Latitude = 16 deg. Longitude = 85 deg 
Altitude = 36km 

Cost functional: {a — atrimY 

This problem has been solved using trial version of General 
Purpose Optimization Software (GPOPS-II). Figures 4 to 
5 show altitude and speed variation. Figure 6 shows the 
ground track, indicating the terminal point reached. 7 
and 8 represent control angle variation to generate this 
trajectory. 

7. PROBLEM 2 DEFINITION 

The reentry vehicle in its journey to ground passes through 
various aerodynamic regimes. Hypersonic phase consid¬ 
ered in last problem is one of them. The aerodynamic 
coefficients of the vehicle vary in these phases. In order to 



Fig. 5. Speed vs Time 



Fig. 6. Ground track 



Fig. 7. Angle of attack vs Time 



Fig. 8. Bank angle vs Time 

simulate and optimize the full trajectory the variation of 
these coefficients should be known and taken into account 
for the optimization problem. In this section the problem 
of reentry of Apollo capsule is considered with following 
properties 

Vehicle properties: 

Mass = 5809 kg, Area = 12.017 sq.m 
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Fig. 9. Cl fit for Mach 15.8 



Fig. 10. Cd fit for Mach 15.8 

8. AERODYNAMIC COEFFICIENTS EOR APOLLO 


Unlike last problem the aerodynamic coefficients are ob¬ 
tained here from the wind tunnel testing report for Apollo 
vehicle . In order to form an analytical expression of the 
variation of coefficients with angle of attack the data from 
report is extracted and the a polynomial curve is fitted 
using Matlab® curve fitting toolbox through the extracted 
datapoints. The entire regime is divided into 5 phases. Eor 
each section it is assumed that the coefficients in that 
particular phase are constant, corresponding to those at 
a particular Mach number in that range. The variation 
of Lift and Drag coefficients is obtained around the trim 
angle of attack in that phase. The equations 29 to 38 show 
such analytical expressions. The figures 9, 10 represent the 
fitted curve against the datapoints for a particular case of 
phase 1, Mach 15.8. 

8.1 Aerodynamic coefficients 

Phase 1 : (Mach >12) Trim angle of attack = 158 degree 

Cl = -0.0042((q; - 167.6 )- 0.033 * ((a... 

-167.6)/9-6)^ - 0.1567((q; - 167.6)/9.6) (29) 

+0.2422 

Cd = 0.0028((a - 164.6)/9.22)3 - 0.035((a... 

-164.6)/9.22)2 + 0.0943((a - 164.6)/ (30) 

9.22) + 1.194 

Phase 2 : (6 < Mach < 12) 

Trim angle of attack = 158.7 degree 


Cl = -0.0027((a - 170.5)/11.52)^ - 0.008 * ((a... 

-170.5)hl-52)^ - 0.017 (q; - 170.5)/... (31) 

11.52)2 - 0.2356((q; - 170.5)/11.52) + 0.1961 

Cd = -0.00068((q; - 160.9)/6.33)^ - 0.001 * ((a... 

-170.5)/11.52)3 - 0.013((q; - 160.9)/ ... (32) 

6.33)2 _ o,1259((q; - 160.9)/6.33) + 1.35 

Phase 3 : (2 < Mach < 6) 

Trim angle of attack = 159.6 degree 

Cl = -0.00041((a - 174.7)/15.5)^ - 0.0255 * ((a... 

-174.7)/15.5)^ - 0.028((a - 174.7)/15.5)2 ... (33) 
-0.3458((a - 174.7)/15.5) + 0.1088 

Cd = -0.0023((q; - 174)/15.58)^ - 0.0018 *((«... 

-174)/15.58)^ - 0.116((q; - 174)/15.58)2 ... (34) 
-0.0867((a - 174)/15.58) + 1.484 

Phase 4 : (1 < Mach < 2) 

Trim angle of attack = 156.3 degree 

Cl = 0.0237((a - 168.5)/14.23)® - 0.0645 * ((a... 

-168.5)/14.23)^ - 0.02355((a - 168.5)/ ... , . 

14.23)^ - 0.06653((a - 168.5) - ... 

/14.23)2 - 0.3706((a - 168.5)/14.23) + 0.251 

Cd = -0.0251((q; - 163.6)/14.74)4 - 0.049 

-163.6)/14.74)3 - 0.003345((q; - 163.6)/... (36) 
14.74)2 - 0.085((a - 163.6)/14.74) + 1.273 

Phase 5 : (0 < Mach < 1) 

Trim angle of attack = 165 degree 

Cd = -0.023((a - 173.8)/15.16)® - 0.0111 ((a... 
-173.8)h5.16)'‘ - 0.07223((a - 173.8)/... 

15.16)3 - 0.02731((a - 173.8)/15.16)2 - ... (37) 

0.00527((a - ... 

173.8)/15.16) + 1.009 

Cd = 0.00227((q; - 164.6)/9.218)3 - 0.035((a - ... 

164.6)/9.218)2 + 0.0943((a - 164.6)/9.218) (38) 

+1.194 

Initial conditions: 

Latitude = —7.67 deg, Longitude = 90.3 deg 
Altitude = 100 km^ Heading = —11.75 deg 

Speed = 8 kmjs., 

Flight path angle = —1.05 deg 
Constraints: 

Angle of attack = —30 deg to 30 deg 
Bank angle = —90 deg to 90 deg 

Linkage constraints connecting each state^ control 
variable in successive phases. 

Terminal condition: 

For phase 1 to 4, the boundary conditions of mach 
number for that particular phase. 

For phase 5, flight Path Angle = 89 deg 

Cost functional: 

Terminal cost : {{Longitude — 85)^ + {Latitude — 16)^) 
Continuous cost : {a — atrimY 
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Fig. 11. Altitude vs time 



Fig. 12. Speed vs time 



Fig. 13. Longitude vs Latitude 

9. PROBLEM 2 RESULTS 


This problem also has been solved using trial version 
of General Purpose Optimization Software (GPOPS-II). 
Eigures 11 to 16 show variation of states. Eigure 17, 18 
represent control angle variation to generate this trajec¬ 
tory. It shows that the capsule reaches close to the desired 
terminal point of 85deg longitude and 16deg latitude. 

10. CONCLUSION AND EUTURE WORK 

Within a limited deviation from reentry initial conditions, 
the terminal point is reached by capsule, while minimizing 
the cost functional. Also the multiphase trajectory opti¬ 
mization is also demonstrated based on the wind tunnel 
data of Apollo vehicle. The current work indirectly ac¬ 
counts for the fuel consumption through deviation from 
trim angle of attack. Also, the consumption due to banking 
manoeuvre is not considered. This is because the current 
optimal control problem formulation is only considering 
translational dynamics. Hence, modelling vehicle attitude 



Eig. 14. Elight path angle vs time 



Eig. 15. Heading vs time 



Eig. 16. Mach number vs time 
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Eig. 17. Angle of attack vs time 

dynamics and reaction control system, and then solving 
the actual fuel minimization problem forms as the poten¬ 
tial future work. 
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Fig. 18. Bank angle vs time 
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